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We investigate a possibility of estimating mass of an isolated rapidly rotating neutron star (NS) 
from a continuous gravitational wave (GW) signal emitted by the NS. When the GW passes through 
the gravitational potential of the NS, the GW takes a slightly longer time to travel to an observer 
than it does in the absence of the NS. Such a time dilation effect holds also for photons and is 
often referred to as the gravitational time delay (or the Shapiro time delay). Correspondingly, the 
phase of the GW from the NS shifts due to the Coulomb type gravitational potential of the NS, 
and the resulting logarithmic phase shift depends on the mass, the spin frequency of the NS, and 
the distance to the NS. We show that the NS mass can, in principle, be obtained by making use 
of the phase shift difference between two modes of the continuous GW such as once and twice spin 
frequency modes induced by a freely precessing NS or a NS containing a pinned superfluid core. We 
estimate the measurement accuracy of the NS mass using Monte Carlo simulations and find that 
the mass of the NS with its ellipticity 10 -6 at 1 kpc is typically measurable with an accuracy 20% 
using Einstein Telescope. 


I. INTRODUCTION 

Mass is a fundamental quantity of an astronomical object. In the case of a neutron star (NS), mass is 
important, among other aspects [l| , as it gives us a clue to study the equation of state of dense matter beyond 
the nuclear density, which is still uncertain through theoretical studies or experiments on the Earth (See, e.g., 
0 )- 

Almost all the measurements of NS masses so far are limited for NSs in binaries. An exception is for the 
nearby isolated NS RX J 185635-754 Q, where its radiation radius and gravitational redshift have been measured 
modulo detailed modeling of its atmosphere 3 (See also the comments mentioned in the Sec. 4 of [5| and Fig. 
3 there). Since an isolated star may have a different evolutionary history from a star in a binary (namely due 
to their mutual mass transfer), the former may follow a different mass function than the latter and hence it 
is interesting to measure masses of as many isolated NSs as possible. This paper proposes a new method to 
measure a mass of an isolated rapidly rotating NS that emits gravitational wave (GW). 

We expect direct detection of GWs within the next decade [6|. The second-generation gravitational wave 
detectors such as KAGRA0, advanced Laser Interferometric Gravitationalwave Observatory (aLIGO) Q, and 
advanced Virgo (aVirgo) [9| are under construction. Conceptual studies of third-generation gravitational wave 
telescopes are also on-going EIU]. One of the most important family of GW sources for those GW detectors 
is a compact binary coalescence (CBC). Detection of GW from a CBC event enables us to measure masses of 
compact stars in the binaries to sub-percent level. A possibility is actively discussed where we can study internal 
structure of the compact stars (if they are regular ones) through the tidal disruption of the binary stars during 
the binary merger ([12|] for a review). However, again, the estimates of masses possible from detection of GWs 
from CBC events are only for stars in binaries. We here focus on another type of GW sources, namely, isolated 
rapidly rotating non-axisymmetric NSs and show that we can estimate masses of isolated NSs by detection of 
those sources with one of proposed third-generation telescopes, namely, Einstein Telescope (ET) ‘la¬ 
in fact, our method proposed in this paper utilizes the fact that phase of GW from an isolated rapidly rotating 
compact star is modulated by the Coulomb type gravitational field of the star itself. This logarithmic phase 
shift is well-known in quantum mechanics fl3l ]. and partly incorporated in the GW waveform from CBC dl- 
This Coulomb type phase shift corresponds to the famous gravitational time delay (or the Shapiro time-delay 
when it is observed in a binary system). But this phase shift seems neglected so far in the context of the 
studies searching for GWs from isolated pulsars. Then our idea to extract mass Mns of an isolated pulsar at 
distance r from GW observation is as follows. In the case of GW, the Coulomb phase shift take a form of 
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^(wgw) = 2w gw MNS ln( 2 w gw r) for GW with frequency w gw [TEj . If we detect two GW modes at frequencies 
Wgw,i and w gw ,2 from the same pulsar, we could find the mass by (5$(wg Wj i) — = 2 oj gWi i AIns l n K 

where K = w g w,i/w g w, 2 - Before examining the feasibility of this idea, we first mention whether a pulsar could 
emit continuous GW at (more than) two frequencies. 

A rotating tri-axial object with its rotational frequency w ro t emits GW mainly at frequency 2w ro t when the 
rotational axis of the object is aligned to the one of the principal axis of its moment of inertia. If the former 
axis is not aligned to any one of the latters and this object rotates as a rigid body as a whole, this object will 
freely precess and emits GW at frequencies close to w ro t and 2w ro t 0]. When the former axis is not aligned 
to any one of the latters and this object has a pinned superfluid layer(s), Jones proposed a possibility that this 
object could steadily rotate at a constant frequency and emits GW at frequencies close to w ro t and 2w ro t PI- 
The difference between the two scenarios 0] and m is that the latter does not freely precess and hence it 
shows no electromagnetic signature of free precession. 

Observational evidences of existence of freely precessing NSs are recently reported by Makishima et al. 0111 
where they found phase modulations in pulsations of the magnetar 4U 0142+61 and IE 1547.0-5408 using the 
Suzaku X-ray observatory. They interpreted that those phase modulations are due to free precessions of those 
objects. There are several other objects that show signatures of free precessions, but it is not clear that current 
theoretical understandings of NS interior structures allow a free precession sustainable long enough so that we 
have any practical chance to observe it at all (See, e.g., fjtil [ 22 } for discussion and references therein). This in 
turn means that the findings by Makishima et al. and further detections of freely precessing NSs give us an 
insight onto the interior structure of NSs and drive our theoretical understandings on it. 

Indeed, majority of pulsars do not show any clear signature of free precession. However, Jones pointed out 
that a NS may have a pinned superfluid. Jones then advocated a theoretical possibility that a NS could rotate 
steadily even when the axis of the moment of inertia of its solid outer crust is misaligned to its rotational 
axis [13. As a result, such a NS emits GW at two frequencies close to w ro t and 2w rot without showing any 
electromagnetic signature of free precession. This proposal by Jones motivates Bejger and Krolak 0] to study 
a search method for GWs from known pulsars at once and twice their spin frequencies. Motivated by these 
observations and the theoretical proposal by Jones, we here assume a pulsar that emits GW at frequencies 
w gWj i ~ w ro t and w gWi 2 = 2 w gWj i ~ 2 w ro t to study feasibility of our idea of estimating the pulsar mass using 
Coulomb phase shifts. The GW radiation at the frequency w gWj i ~ lo tot is sometimes called a “wobble mode” or 
“wobble radiation” for a freely precessing NS. Here we would like to include a non-precessing NS that emits GW 
at that frequency, we here call this mode “the first harmonic mode” while the GW at 2 oj “the second harmonic 
mode” (or “the first overtone mode”). Applications of our idea to GW modes other than the ones considered 
here or extensions to the case where multiple GW modes more than two are available should be straightforward. 

The rest of this paper is organized as follows. In Sec. El we construct a GW subject to a gravitational 
Coulomb phase shift due to the gravitational potential created by the GW source. In Sec lIIIl we apply our 
formulation to an isolated NS and examine measurement accuracies of the NS parameters. Summary and 
conclusion are presented in Sec lIVI Throughout this paper, we adopt the geometric unit system in which both 
the light speed c and the gravitational constant G are unity. 


II. GW WAVEFORMS PROPAGATING THROUGH NS GRAVITATIONAL POTENTIAL 

A. Green’s function 


In this section, we derive the explicit form of the gravitational phase shift for GW caused by the gravitational 
potential of the source itself. The derivation in this section is basically based on the section 2 in [15|. Consider 
an isolated GW source with mass M and located at r = 0. The dominant static part of the exterior gravitational 
held far away from this object is well approximated by 


ds 2 = - 1 - 


2 M 




dt 2 + ^1 H-^ 5ijdx l dx-i. 


(1) 


A metric perturbation h ,iU = r 7 ^ — sj—gg^ with respect to this metric satisfies the following wave equations 

4 M 


- 1 + 


d + 


h^ = 0 


( 2 ) 


outside of the GW source in the harmonic gauge h ,iv = 0 where rf* 1 ' denotes the hat spacetime metric. 
Motivated by this observation and to capture the dominant effect of the gravitational held of the source onto 
the GW phase, we express the full Einstein equations including the source under the harmonic gauge condition 
in the following form, 


□ M h^ v = - 16 ? TT^. 


( 3 ) 
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The effective energy-momentum pseudotensor rand the differential operator Dm are defined as 


□m = 
_ 

_ 

•-M — 


- 1 + 


4 M\ 


d 


r / 

+ C LI 

h^ a 8 h vf) a + 2WhP v o j + b? k h? v ik + 


{-g) (T^ + + O , 


1 

107T 




(4a) 

(4b) 

(4c) 


where denotes the Landau-Lifshitz pseudotensor [0 and M denotes the total mass of the source. It should 
be noted thcit the term nppenrs in the both sides in Ep. ©■ Because the 00 component of 

the metric perturbation is evaluated in the wave-zone as h 00 = 4M/r to the lowest order in the weak field 
approximation, we intentionally added this term to the both sides in Eq. ©■ Introducing the “potential 
barrier” 4 M/r in the wave operator, we can take into account the deviation of the curved spacetime light cone 
from the flat spacetime one to the lowest order in the weak field approximation. As will be seen in the end of 
this section, this mass term brings the gravitational phase shift in the GW phase. 

Equation. © can be formally solved by using the retarded Green’s function G^ (x,x ') as, 


h M!7 (x) = —167T 


d 4 x 7 G^ (x, x') (x 7 ), 


(5) 


where the retarded Green’s function obeys the following equations. 

□mGm +) (x, x') = <5 (4) (x - x'), (6a) 

G<+ } (x, x') = J2 [ dco e^sgn (w) [<F + ^ m (x) ( x ') H (r - r') + F Suiem (x) ( x 7 ) H (r 7 - r)] . 

tm 

(6b) 


The functions and 4' SWm 

0 : 


are composed of the spherical Coulomb wave functions ( p ) and F. ( p ) 


y SbJ * m (x) 



-iut -1 (+) 

P u l 


( 9 ) Y tm ( 6 », <j >), 



P~ X F t (p) Y tm ( 9 ,0) 


(7a) 

(7b) 


where ae = argr (£ + 1 — 2 iMui), p = rw, and H (x) denotes the Heaviside step function. The radial functions 
u i + ' > (p) an d Ff„ ( p ) can be evaluated in the wave zone or the near zone as follows. 


4 +) (P) 

Mp) 


exp 


i [ p + 2Mu>hi2p — — £tt 


for 


cep 


i+i 


for p —> 0, 


(8a) 

(8b) 


with ce = 2 t ’e nMu |T (£ + 1 — 2iMco)\ / (2 i + 1)!. Since the distance to the source is much greater than the size 
of the source r r 7 , substituting Eas. (17al) - (18bl) into Ea. (l6bl) brings the Green’s function to the following form, 


G[, +) (x,x 7 )=£ 
i 


(~i) e n (t> x' (L) 
2 ( 2 tt F.fr 


dcoT{l+£- 2iMu) e *M Uuj e e -i(t-r-2Mln2 U r-t') U + Q , 


(9) 


where the angle bracket attached to the indices denotes the symmetric-trace-free operation, the suffix L stands 
for l tensorial products of vectors, and the unit vector n? is defined as = x J /r. Combining Eqs. © and (FUl . 
we obtain the GW waveforms affected by the gravitational potential of the source. The authors of Ref. [15| 
then adopted slow motion approximation and expanded this Green’s function in Mto and arrived at the GW 
waveform including the tail term which was first derived by Blanchet and Damour fl4| . They then pointed out 
that the dominant part of the tail term originates from the deviation of the true light cone from the flat light cone 
due to the gravitational potential of the GW source itself, which appears as the phase shift 2iMoj\n2u}r + iae 
in Eq. ©• Since we focus on this gravitational phase shift for GWs, we drop the Moj term in the amplitude 
but retain the Mw term in the phase up to the first order in our analysis. 

Finally we note that this phase shift corresponds to the well-known gravitational time delay (or the Shapiro 
time delay in the case of a binary). The GW propagates along the true light cone of the curved spacetime, and 
so takes a slightly longer time to travel to the detector than it does in the absence of the source. 
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B. Continuous GWs through NS gravitational potential 

Let us now assume an isolated rapidly rotating NS with mass Mns as a GW source. As mentioned in the 
introduction, an isolated rapidly rotating NS would emit GWs at distinct frequencies at the same time. We 
define the volume integral of the stress-energy tensor as, 

S jk (t) = j d 3 x T jk (t, x) . (10) 

If the source has different frequency modes oj\, ui 2 , • • •, the Fourier component of S^ k ( t ) can be decomposed into 
the sum of the delta function, 


S* ( W )=X> sifS^-Wn), 

n 


(ii) 


where is the coefficient corresponding to the 71 -th mode. Combining Eqs. ©, ©-(ED- we arrive at the 
GW waveform including the gravitational phase shift in the quadrupole approximation, 


h? k (t) = j J2^n k e~^ nit \ 

( 12 a) 

(t) = it - r - 2 M ns [ln( 2 rw n ) + C]) u n . 

( 12 b) 


The constant C is the Euler’s number. 

It is easy to see that when Mns 0, Eq. (I12al) with Eq. (I12bl) reduces to the well-known quadrupole grav¬ 
itational waveform as expected. Inclusions of higher order mass multipoles or current multipoles, if necessary, 
are straightforward. 


III. MEASUREMENT ACCURACY OF NS MASS 

In this section, we apply the GW waveform including the gravitational phase shift effect to an isolated rapidly 
rotating NS that emits GWs at two frequencies, w gWi i = u and w gWi2 = 2ui. 


A. GWs from an isolated rapidly rotating non-axisymmetric NS 


A non-axisymmetric NS rotating around its principal axis produces the quadrupole GW radiation with a 
frequency 2w rot , where w ro t denotes the rotational frequency of the NS. However, the angular momentum axis 
of the NS does not gen erally align with its principal axis because of, for example, the strong toroidal magnetic 
held inside the NS [25(. When the NS deviates from the spherical symmetry, the misalignment between the 
angular momentum axis and the principal axis in general leads to a free precession of the NS. It would also 
be possible that the NS steadily rotates at w ro t at the fixed rotational axis in an inertial frame if it has a 
pinned superfluid component inside it G2- Those two kinds of NSs radiate the GW with two frequency modes 
w™ 1 = uj and w gWi2 = 2 ui which we call the first harmonic mode and the second harmonic mode, respectively 


The GW waveforms in these two scenarios are generally different from each other for tri-axial NSs |16|, [17], 
!H, [ 2 ^, [ 23 . ■ Interestingly, in the case of bi-axial NSs, those two families of the GW waveforms take the same 
form (shown e.g., in [28(), which we adopt in this paper to apply our result to both scenarios. 

We then adopt the following GW waveforms 


h + = A +i 1 Cos'Ll + A +: 2 cos’ll, 

(13a) 

h x = A Xi i sin 'F 1 + A Xj2 sin'll 

where A and T denote the amplitude and the phase and are defined as 

(13b) 

A + .i = —ho sin 29 sin l cost, 

(14a) 

A +i 2 = sin 2 8 (l + cos 2 t) , 

(14b) 

A x .i = -/iosin20sini, 

(14c) 

A x ,2 = ho sin 2 0 cos 1 , 

(14d) 

’Ll = (t — r — 2Mns In 2 rut) lo + ipo > 

(14e) 
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'$ r 2 — 2 (t — r — 2Mns ln4rw) ui + 2ipo- (14f) 

The angle i is called the inclination defined as the angle between the rotational axis and the line-of-sight, and 
9 is the misalignment angle between the principal axis and the angular momentum axis. The angle ip® is the 
GW phase at the reference time t = 0 at the NS [32j. The overall amplitude ho is given by ho = Aeluj 2 /r where 
I and e are the principle moment of inertia and the ellipticity, respectively. Note that the log-term appears in 
the GW phase due to the gravitational potential of the NS. 


B. How to extract the NS mass from GW signals 

In this paper, we take an example of GWs from an isolated rapidly rotating non-axisymmetric NS to demon¬ 
strate the measurement accuracy of the NS mass. A NS would emit GWs of several modes such as the first and 
the second harmonic modes mentioned in this paper, wobble mode if the spin axis of the NS precesses, /-mode 
and r-mode if the corresponding stellar oscillations are excited. If the initial GW phases of two modes at the NS 
are known, or the ratio of the two GW phases is equal to the ratio of the GW frequencies of the two mode, the 
NS mass can be estimated by subtracting the GW phase of one mode from the other with the ratio multiplied 
appropriately. In the case of the first and second harmonic modes considered in this paper, this assumption is 
guaranteed because the both modes are generated by the deformation of the crust of the NS. Multiplying 
by 2 Subtracting Eq. m from Eq. (THcl) . we obtain 

2'f'i - = 4 M ns w In 2. (15) 

Equation. (fl5l) indicates that the difference between the two modes separates the NS mass and the distance 
to the NS. Therefore, the NS mass alone can be estimated by GW observations. However, the accuracy of the 
constant phase is typically of the order of 0.1 for a signal to noise ratio (SNR) 10 even in targeted searches for 
known pulsars where w ro t(— w gWj i) is well-known by electromagnetic observations. Then, the second generation 
GW detectors would have a little chance of detecting the log-term appeared in Eq. (fl5l) with enough sensitivity 
to estimate the NS mass. On the other hand, a third generation GW detector whose sensitivity is about an 
order of magnitude better than that of the second generation ones could determine the NS mass from the GW 
signal. 

As will be discussed in the next section, we consider targeted searches for GWs from “known” rapidly 
rotating non-axisymmetric NSs by a third generation detector. Those objects are “known” in the sense that 
electromagnetic observations and/or the second generation detectors will have detected them before the search 
using third generation detectors. We then assume that the positions of the NSs in the sky and its GW frequencies 
(or rotation frequncies) and their times derivatives are known in advance. Those parameters are not search 
parameters in the following hypothetical search with a third generation GW detector. 


C. Monte Carlo simulations 

In this section, we investigate how accurately the mass of the NS can be determined by GW observations. 
Measurement errors are estimated by using the Fisher analysis which is briefly summarized in Appendix [Bj A 
GW signal from an isolated rapidly rotating non-axisymmetric NS is given by Eqs. (lAHl-dAfil). The GW signal 
is characterized by 6 waveform parameters A = {ho, 9, t, ip, epo, -A/ns} in targeted searches. To demonstrate the 
accuracies of the waveform parameters A, we perform Monte Carlo simulations for the observation time T 0 b s = 3 
years in which we treat a, <5, 9, l, and ip as random variables. Here, the sky position of the NS is specified by 
two parameters, the right ascension a and the declination <5. The NS mass Mns is assumed to be 1.4 M 0 . We 
adopted e = 10 -6 [2j|, I = 10 38 kg • m 2 , and selected distances and frequencies to compute the amplitude ho- 
For each simulation, we generate 10000 sets of random variables each of which is distributed according to the 
uniform distribution. 

Figure. [T] displays the cumulative distribution functions of relative errors in the NS mass estimation 
AMns/AIns from our Monte Carlo simulations for ET observations. We investigated the cumulative distri¬ 
bution functions for two different GW frequencies 300 Hz and 500 Hz. These figures show, for example, in the 
case of the GW signal with / rot ~ / gw ,i = w gWj i/(27r) = 500 Hz from NSs at r = 1 kpc in ET observations, 
that masses of about a 70 percent of isolated rapidly rotating non-axisymmetric NSs are measurable with an 
accuracy of AMns/A/ns — 0.2. It can be detected with SNR = 0(10) by a single second generation detec¬ 
tor because the sensitivity of ET is about an order of magnitude better than that of second generation one. 
Therefore, if an isolated rapidly rotating non-axisymmetric NS is detected by a network of second generation 
detectors, the mass of the NS can be determined with an accuracy of AMns/ATns = O (0.1) in a single third 
generation detector. Also, Fig[l] indicates that the accuracy of the NS mass becomes better and better as the 
GW frequency increases. This behavior can be traced to Eq. m■ The effect of the gravitational phase shift 
is proportional to the rotational frequency. So the NS mass can be more easily extracted from the GW signal 
with the higher frequency. 
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Figure. [2] shows that the scatter-plot of the NS mass accuracy AA/ns/^ns as a function of the inclination l 
(left) and the angle 9 (right) for GW signals with / gw ,i = 500 Hz from NSs at r = 10 kpc. The left panel of 
Fig. [2] indicates that the NS mass cannot be estimated at all at i = 0,7r while well estimated at around t = 7r/2. 
This is because the amplitudes of the first harmonic mode in both the plus and cross polarizations vanish at 
i = 0, 7 r from Eqs. (I14al) and (I14cl) . In these cases, since the GW signal has only the second harmonic mode, 
the effect of the gravitational phase shift is unobservable. The right panel in Fig. [2] indicates that the NS mass 
cannot be determined for 0 = 0,7r/2. This fact can be interpreted in the same way as the left panel of Fig. [2] 


(a) f = 300 [Hz] 



(b) f = 500 [Hz] 



Relative error of mass AM/M 


Relative error of mass AM/M 


FIG. 1: The cumulative distribution function as a function of relative error of NS mass AMns/AIns for two different 
GW frequencies (a) 300 Hz and (b) 500 Hz. The solid line, the dashed line, and the dashed-dotted line corresponds to 
NSs at the distances r = 1 kpc, 10 kpc, and 50 kpc respectively. For NSs at the distance of r = 50 kpc that emit GW 
at the first harmonic GW mode frequency of / g w ,i = 300 Hz, measurement accuracy becomes more than 100%. 




FIG. 2: The scatter-plot of relative error of NS mass AAIns/A/ns an d inclination i (left), misalignment angle 0 (right). 
The GW signals are assumed to have the GW frequency 500 Hz from NSs at r = 10 kpc. 


IV. CONCLUSION 

In this paper, we have focused on a gravitational phase shift for a GW from an isolated rapidly rotating 
non-axisymmetric NS. When the GW passes through the gravitational potential created by the NS, the GW 
phase is stretched due to the gravitational phase shift effect depending on both its rotational frequency and its 
mass. We have constructed the explicit form of the resulting phase shift based on Ref. [u| and obtained Eq. 
(I12bl) . We have shown that if the GW from the NS has different frequency modes such as the first and the 
second harmonic modes as proposed in the mass of the isolated NS can be, in principle, determined 

by making use of the phase difference of the different modes. We have given an example of an isolated rapidly 
rotating non-axisymmetric NS to estimate the measurement accuracy of the NS mass using Fisher analysis in a 
targeted search. Performing Monte Cairo simulations in 3-year ET observation, we have obtained a cumulative 
distribution function of a relative error of the NS mass AMns/ATns for different GW frequencies 300 Hz and 
500 Hz and different distances of r = 1,10, and 50 kpc as shown in Fig. |1] The range of the NS mass can be 
typically delimited with an accuracy of AMns/ATns = O (0.1) with third generation detectors for such a GW 
signal that can be detected by a second generation GW detector with SNR = 0(10). 
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Appendix A: GW signals 


The GW waveforms from an isolated rapidly rotating non-axisymmetric NS are presented in Sec. IIII Al As 
can be seen in Eq. (U3al) and (113bl) . the GW signals consist of two components with frequency / gw ,i and 2/ gWi i. 
It is convenient for our purpose to decompose the GW signals as follows [28| . 


The amplitudes Ajk are defined as 
An = ho sin 28 
A 12 = hg sin 29 
A 13 = h 0 sin 28 
A 14 = ho sin 29 


2 4 

h{t) =J2J2 A i kh i k (*)' 

j=i fc=l 


— sin 2 l cos 2'i p cos 0 o-sin l sin 20 sin 0 o 

8 4 

— sin l cos 20 sin 0 o + — sin 2 t sin 20 cos 0 o 

4 8 

— — sin 2 l cos 20 sin 0 g — — sin l sin 20 cos <j 


— sin l cos 20 cos 0 o — — sin 2 l sin 20 sin <; 


A 2 i = ho sin 2 8 
A 22 = ho sin 2 8 
A 23 = ho sin 2 8 
A 2 4 = ho sin 2 8 


(l + cos 2 i) cos 20 cos ( 2 </>o + 0 S ) — cos l sin 20 sin ( 20 O + 0 S ) 
(1 + cos 2 l) sin 20 cos ( 20 o + 0 s) + cos t cos 20 sin ( 20 o + 0 S ) 


— - (l + cos 2 t) cos 20 sin ( 20 o + 0 S ) — cost sin 20 cos ( 20 o + 0 S ) 

— i (l + cos 2 t) sin 20 sin ( 20 o + 0 S ) + cos t cos 20 cos ( 20 o + 0 S ) 


(Al) 

(A2a) 

(A2b) 

(A2c) 

(A2d) 

(A2e) 

(A2f) 

(A2g) 

(A2h) 


where ho,i, 0,0 and 0 o are the overall amplitude, the inclination, the polarization phase, and the initial phase, 
respectively. The phase 0 S expresses the gravitational phase shift due to the gravitational potential of the NS 
and corresponds to Eq. IT5l) . It is described by 

0s = 87rM NS /gw,i In 2. (A3) 

The time-dependent parts hjk (t) in Ea. dAll) characterize the shape of the GW signal and are defined as 


hji (t) = a (t) cos [ 7 $ (i)], 
h j2 (t) = b(t) cos \j$(t )], 
hj 3 (t) = a (t) sin [jd> (t)], 
hji (t) = b(t) sin [j$(t)] 


(A4a) 

(A4b) 

(A4c) 

(A4d) 


where the function a (t) and b (t) are modulation amplitudes given in [28( . The GW phase after the barycentric 
corrections is expressed by 


$ (t) = 2nf gwA (t + n-r (t)) (A5) 

where n denotes a unit vector pointing toward the NS from the solar system barycenter (SSB), and r (t) denotes 
the separation vector pointing from the detector to the SSB. Note that we neglected the spin-down effect in the 
GW frequency for simplicity. The decomposition of the GW signals in Eq. m separates known parameters 
and unknown parameters in targeted searches. 






















Appendix B: Fisher analysis 


In this appendix, we give a brief review of the Fisher analysis |30j, 31]. The detector output s (t) can be 
expressed by linear sum of the GW signal h (t) and the detector noise n(t), s ( t ) = h{t)+n ( t ). When the noise 
is stationary and obeys the Gaussian distribution, measurement errors of the waveform parameters A are found 

by 

("T„,= V« A W> = v'd'-'t, (bi) 


for large signal-to-noise ratio where (AA J ) rms is the root-mean square of the waveform parameters AA J . The 
matrix Tis called the Fisher information matrix and is defined by 


r jk — 



(B2) 


where (-|-) denotes the noise-weighted inner product, and djh denotes the derivative of the GW signal with 
respect to the waveform parameter A j. For instance, the inner product between two time functions x(t) and 
y (t) is expressed by 


(x\y) = 4Re 



x{f) V* (/) 
Sn(f) 


(B3) 


where S n (f) is the one-sid ed p owerspectral density of the detector. In this paper, we have adopted S n (f ) of 
ET given in the Table 1 of (32|. 

When the frequency of the signal is nearly constant, it is convenient to introduce a new inner product 

2 pT, obs/2 

01 \y) = 7 ^— x (t) y (t) dt (B4) 

Jobs J-T obs /2 

where T 0 b s denotes the observation time [28|. Since the time scales of change in modulation amplitudes a (t) and 
b (t) are of the order of the time scale of the Earth’s rotation, the GW phase $ (t) oscillates much more rapidly 
than modulation amplitudes. Then, the inner product between the amplitudes hj k can be well-approximated 
as 


- (Ml h os) - l Hl a ) > 

(B5a) 

(MIM - (MilM - \ ( b W b ) i 

(B5b) 

(MUM) - (MIIM) - \ ( a W b )- 

(B5c) 


Combining Eqs. (IB2l) - (IB5cl) . we obtain the Fisher matrix in our analysis, 

T 0 bs 


M = E 


(djAadkAn + djA t3 d k A a ) (a||a) + ( d j A l2 dkAi 2 + djAudkA e4 ) (b\\b) 


=1 2 ^‘ 

+ {djAnd k A£2 + djA^dkAn + djA^dkA^ + dj Ap 4 dkA^) 00 ) 


(B6) 
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